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ABSTRACT 

We apply the empirical method built for redshift z = in the previous work of 
Wang et al. to a higher redshift, to link galaxy stellar mass directly with its hosting 
dark matter halo mass at redshift of around 0.8. The M stars -M in f a u relation of the 
galaxy stellar mass M star s and the host halo mass Mi n f a u is constrained by fitting 
both the stellar mass function and the correlation functions at different stellar mass 
intervals of the VVDS observation, where M in f a u is the mass of the hosting halo 
at the time when the galaxy was last the central galaxy. We find that for low mass 
haloes, their residing central galaxies are less massive at high redshift than those at 
low redshift. For high mass haloes, central galaxies in these haloes at high redshift are 
a bit more massive than the galaxies at low redshift. Satellite galaxies are less massive 
at earlier times, for any given mass of hosting haloes. Fitting both the SDSS and 
VVDS observations simultaneously, we also propose a unified model of the M stars - 
Mi n f a u relation, which describes the evolution of central galaxy mass as a function of 
time. The stellar mass of a satellite galaxy is determined by the same M sta rs-Mi n f a ii 
relation of central galaxies at the time when the galaxy is accreted and becomes a 
sub-component of a larger group. With these models, we study the amount of galaxy 
stellar mass increased from z^ 0.8 to the present day through galaxy mergers and star 
formation. Low mass galaxies (< 3 x 10 w h~ 1 M@) gain their stellar masses from z^ 0.8 
to z — mainly through star formation. For galaxies of higher mass, we find that the 
increase of stellar mass solely through mergers from z = 0.8 can make the massive 
galaxies a factor ~ 2 larger than observed at z = 0, unless the satellite stellar mass is 
scattered to intra-cluster stars by gravitational tidal stripping or to the extended halo 
around the central galaxy that is not counted in the local observation. We can also 
predict stellar mass functions of redshifts up to z ~ 3, and the results are consistent 
with the latest observations. Future more precise observational data will allow us to 
better constrain our model. 

Key words: galaxies: masses - galaxies: high-redshift - galaxies: haloes - cosmology: 
dark matter - cosmology: large-scale structure 



1 INTRODUCTION 

To study how galaxies form and evolve in their hosting 
dark matter haloes, a lot of efforts have been made to 
link galaxy properties with the properties of dark matter 
haloes which they reside in. The usual meth ods used in- 
clude galaxy kinematics llErickson et all Il987t) and galaxy 
lensing( M andelbaum et al .1 120051 , 120061 ). which measure the 
mass of hosting dark matter haloes directly. Semi-analytic 
models trace the gas cooling, star formation, and feed- 
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back processes 'ab initio' to get the properties of galaxie s 
of present dav (|de Lucia fc Blaizoll2007l : iBower et aij|2006h . 
Halo occupation distribution models study the galaxy- 
halo connection empirically, to model galaxy properties 
using certain assumed formu l a to describe the galaxy- 
halo relati ondJing et all 1 19981; iBerlind fc Weinberg! | 2002| : 
Yang et al.l|2003l; Ivale fc Ostrikerll2004l ; IConrov etallbood ; 



Wang et al.ll2006t ) 



For the models that describe the formation and 
evolution history of galaxies, the statistics that are 
commonly used to constrain these models include: 
number statistics such as number density, luminosity 
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function and stellar mass functionfeullock et all |2002|; 
IZehavi et all 120051 ; iMoster et all 12003 ) , spatial clustering 
prop er ties described by correlation f uncti ons ljjing et al.l 
1 19981 ; lYang et all |2003|; IZehavi et all 120051) void prob- 
ability distributio nllVale fc Ostrikerl |2004 ). pairwise 
velocity dispersiondjing et all 1 19981 ) . the Tully-Fisher 
relation (|Yang et all 
galaxies and the 
colour ! 



2003), 



Zheng 



the colour distribution of 
cluste ri ng d e pendence on galax; 



P_ 

van den Bosch et all 120031 ; iKravtsov et all 
2004l ; IWang et al.ll2007l ).~ 



alaxy 
2004; 



The current models of galaxy-halo connection mainly 
focus on low redshift study, particularly of the present day, 
simply because we can handle well the observational statis- 
tics only at low redshifts. With the development of large 
scale galaxy redshift surveys, observations are obtained not 
only for the local Universe, but also toward higher red- 
shifts. Surveys aiming at studying the properties of hig h 
redshift galaxies includ e DEEP2 survey l|Davis et all [2003). 
the COMBO-17 surv ey ll Wolf et al.1120041). the VIMOS-VLT 
Deep Survev(VVDS) (lLe Fevre et alll2005h. the Cosmic Evo- 
lution Survey (COSMOS) llScoville et al.N2007l ) and zCOS- 
MOS survev l|Lillv et all l2009h . Withthe help of these 
surveys, luminosity functions of differ ent types of galax- 
ies are obtained up to redshift z > 7 (|Reddv et all 120081 ; 
iBouwens et all 120081) . Correla tion functions in luminosity 
bins reaches redshift of z ~ l jCoil et al. 2006; iPollo et all 
2006). Stellar mass functi ons have been detected for galax- 
ies up to redshift of z ~ 5 l|Drorv et al.ll2005l ; iFontana et all 
120061 ; lElsner et all I2OO8I ). Correlation functions in bins of 
st ellar mass have been stud ied for galaxies of redshift z ~ 
l l|Meneux et al.ll2008l . I2OO9I) for VVDS and zCOSMOS ob- 
servations. 

Based on the observational data obtained at high red- 
shifts, several works have been done to model the proper- 
ties of galaxies at early epoch. Some work s use the HOD 
models to study high z galaxy prop ertiesdCooravl 120051) , 
as well as galaxy clustering properties ( Bullock et al.ll2002l ; 
Yan et al.ll200l; ICoorav fc OuchilE^oll IConrov et ai]|2006l ; 



White et al.ll2007l ). which focus mainly on the c lustering de- 



ende nce of galaxy luminosity. Most recently, IZheng et all 
2007) uses the HOD method to model the clustering of 
DEEP2 galaxies as a function of luminosity, which reaches 
redshift of z ~ 1. While luminosity is the most studied 
and easily got property of a galaxy, stellar mass is never- 
theless a more fundamental property, since luminosi t y may 
be affected a lot by dust attenuation. IMoster et all (2009) 
uses a statistical approach to determine the relation between 
galaxy stellar mass and its hosting dark matter halo, and 
constrains the evolution of galaxy stellar mass by fitting 
the stellar mass fu nctions at different redshifts taken from 
iDrorv et all (|2004l ). 

In the previous work o f lWang et all (|2006h , an empirical 
method has been used to link galaxy stellar mass directly 
with its hosting dark matter halo mass. The method falls in 
between the semi-analytic approach and the halo occupation 
distribution approach. Positions of galaxies are predicted by 
following the merging histories of halos and t he trajectories 
of su bhaloes in the Millennium Simulation (|Springel et all 
120051 ). The stellar mass of galaxies at redshift is related 
to the quantity Mi n f a ii by a double power law function. 
Mi n f a u is defined as the mass of the halo at the time when 
the galaxy was last the central dominant object. Parameters 



describing the function are constrained by fitting both the 
stellar mass function and the correlation function s at dif- 
feren t stellar mass intervals from SDSS observation l|Li et all 
2006). The derived M s tars-M in f a u relation is in excellent 
agreement with the dete rmination from galaxy-gal axy weak 
lensing measurement of Mandclba um et all 1 20061 ). In this 
study, we will apply this method to an earlier epoch. By fit- 
ting the statistic results of VVDS observation, we can study 
the connection between galaxy mass and its hosting halo 
mass at redshift of around 0.8. Based on the relations ob- 
tained both of today and of higher redshift, we can study 
the evolution of galaxy stellar mass from z ~ 0.8 to present 
day. 

After a galaxy falls into a larger group and becomes a 
satellite, its surrounding gas is shock-heated and star forma- 
tion ceases in a short timescale. The stellar mass of satellite 
galaxies should remain about the same as the time when 
they are accreted. In this case, the stellar mass of a satellite 
galaxy is determined by the M a tars-Mi n f a ii relation of cen- 
tral galaxies at the time of infall. This inspires us to describe 
the M s tar s-Mi n f a u relation for all galaxies at any redshift in 
a uniform way, by modelling the evolution of M s tars-Mi n f a u 
relation of central galaxies. Assuming that satellite stellar 
mass does not change after infall, the relation for satellite 
galaxies follows the relation of central galaxies at an ear- 
lier epoch when it is accreted. We will explore this unified 
model in §4. With this model, we can also test if a significant 
amount of satellite disruption by tidal forces is required by 
current observations. 

This paper is organised as follows: in Sec.[5J we present 
the model for fitting VVDS observations to get the rela- 
tion between galaxy stellar mass and the hosting halo mass 
at z ~ 0.8. Based on the models describing galaxy stellar 
masses both at low and high redshifts, we analyse in Sec. [3] 
how galaxies gain their masses from redshift of 0.8 to to- 
day. In sec. 3] we build a unified model to fit observations of 
both low and high redshifts simultaneously, assuming that 
satellite stellar mass is determined by the M s tars-Mi n j a u 
relation of central galaxies at the time of its accretion, and 
study the mass growth of galaxies from z = 0.8 based on 
this model. In sec. [5] we predict the stellar mass functions 
of higher redshifts based on our two best-fit models, and 
compare our results with recent obs ervations. This work is 
based on the Millennium Simulation l|Springel et alll2005l ). 



2 MODEL 



As me ntioned in Sec.[TJ in the previous work of IWang et all 
(2006), the relation between the galaxy stellar mass and its 
hosting halo mass at infall time have been studied by fitting 
both the stellar mass function and correlation functions at 
different stellar mass bins from the SDSS observation. This 
relation can be described by a double power law form for- 
mula. To study the M a tars-Mi n j a u relation at higher red- 
shifts, as a first test, we assume that the M s tars-Mi n f a u 
relation at higher redshifts are the same as that of present 
day. We test whether the resulting stellar mass functions and 
correlation functions are consistent with the observation at 
higher redshifts. W e simply appl y the b est-fit M s t ars -Mi n f a u 
relation at z=0 of Wang et al. (2006) to higher redshifts, 
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Figure 1. Correlation functions derived by applying the Mst ars-Mj n f a u relation of z=0l lWang et al1l2006l) to z=0.83, compared with 
the observational results from VVDS of redshift of around 0. 8 l|Meneux et al]|2008h . Symbols with error bars are from observation, while 
dotted lines represent the model prediction. 



and derive stellar masses of all galaxies at each time. The 
M star s-M i n f a u relation is described as follow: 

2 



M a tar 



M 



-xk. 



The scatter in \og(M s tars) at a given value of Mi„f a ii was 
described with a Gaussian function of a width a. The best- 
fit model to the SDSS observation had the following pa- 
rameters: M = 4.0 x IO'^-'Mq, a = 0.29, = 2.42, 
log k = 10.35 and a — 0.203 for central galaxies and 
M = 4.32 x 10 11 /i" 1 M Q , a = 0.232, f3 = 2.49, logfc = 10.24 
and a = 0.291 for satellite galaxies. 

Once we get stellar masses of galaxies at a certain red- 
shift, we can calculate the stellar mass function and also 
the clustering results at different stellar mass bins at that 
redshift, combined with the dynamical information of sub- 
structures in the simulation. The position of each galaxy is 
derived directly from Millennium Simulation by following 
the evolution of substructures. Fig. Q] shows the projected 
correlation functions at three different stellar mass bins at 
redshift of around 0.8. S ymbols with error b ars are results 
from VVDS observation (|Meneux et al"1l2008l ). and dashed 
lines are the derived results from our simple model. The 
projected correlation functions predicted by the model are 
in reasonably good agreement with the observation from 
VVDS. 

Fig. [5] shows the observed stellar mass functions 
at different redshifts(symbols) , compared with the stel- 
lar mass functions derived from our model(lines) , with 
the M star s-M i n f a u relation taken the same as that of 
present day. Considering the fact that different initial mass 
functions(IMF) were adopted in different observations, we 
conv ert all stellar masses to the case of the Chabrier 
IMF (|Chabrierll2003l 1. in this figure and all the figures of 
stellar mass functions h ereafter. The g alaxy mass obtained 
with the Salpeter IMF rtSalpeterl \l 955) is divided by 1.70, 
and t hat with the KroupajKr oupa 2001) IMF is divided by 
1.104 l|Cowie fc BargerlbOOct ). It is clear from Fig.[2]that the 
stellar mass functions are not reproduced under the sim- 
ple assumption that M 3 tars-Mi n f a u relation does not evolve 
with time. Observationally, the number of galaxies with low 
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Figure 2. Galaxy stellar mass functions at different red- 
shifts. Symbols with error bars are observatio nal results . Blac k 
diomends are SDSS observation of z ~ O.l dLi fc White! |2009|) . 
Black points are VVD S results in the redshift bin of [0.7, 0.9] from 
Pozzetti et a.1.1 J2007f). G reen, red and blue points are results from 



MarchesinietalJ l^ih. in three redshift bins. Lines are model 



prediction, with black, green, red and blue ones corresponding to 
results at redshifts of 0.83, 1.5, 2.07 and 3.06, respective ly. Stellar 
mass es of galaxies are normalized to the Chabrier IMF dChabrierl 
l2003t) . 



stellar mass decreases dramatically towards higher redshifts, 
while the number of high mass galaxies stays roughly un- 
changed with different redshifts. The model derived results, 
however, show an opposite trend. The number of low mass 
galaxies evolves quite slowly, and remains almost the same 
till redshift 2, while the number of high mass galaxies evolves 
a lot, with a much smaller number of galaxies existing at 
higher redshifts. These results show that the M a tars-Mi n f a ii 
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Figure 3. The best-fit model results for fitting both stellar mass function ljPozzetti et al ■120071 ) and correlation functions in three different 
stellar mass binslMe neux et alj .2008) from VVDS observation. Symbols are observational results, and dashed lines are the best fit model 
results. 



relation must vary at different redshifts. A reasonable model 
should in general give more massive galaxies and fewer low 
mass galaxies at high er redshifts than at the local universe. 
As shown in Fig.3 of iMoster et al.l (|2009l ). the stellar mass 
function is more sensitive to the change of model parameters, 
while the \ 2 °f correlation functions stays natter around 
minimum in a large range of parameter sets. This explains 
why the correlation functions can be well reproduced while 
the stellar mass functions show such a large discrepancy. 
Therefore, correlation function alone is not enough to con- 
strain the M a t ar a -Mi n f a ii relation at z ~ 0.8. At the end of 
this section, we will show that stellar mass function alone is 
not enough either to give a good constraint on the relation, 
at least in fitting the current observational data results. 



From Fig. [T] and Fig. [2] we know that we need to fit 
both the stellar mass function and the two point correla- 
tion functions to constrain the underlying M s tars-Mi n f a u 
relation at higher redshifts. Current studies of VVDS ob- 
serva tion give both the stellar mass function l|Pozzetti et al.l 
120071) and correlation f unctions in different stellar mass 
bins (|Meneux et all l200j ) at redshift of ~ 0.8. We there- 
fore focus on constraining Mstars-Mi n f a ii relation by these 
VVDS results, building models based on the si mulation out- 
put at redshift of 0.83. Following the method o f lWang et al.l 
(2006), we assume that the M s t ar s-Mi n faii relation at red- 
shift of ~ 0.8 can be described by a double power law form. 
The relation is determined by five parameters for central 
and satellite galaxies separately, which includes in total 10 
parameters in the modelling. When applying the modelling 
to VVDS results, we notice that compared with the obser- 
vations of the local Universe, the observational results from 
higher redshift survey like VVDS provide much fewer data 
points. Besides, the error bars of these data points are still 
large, which yields weak constraint o n the model. There fore, 
on the basis of our previous model of lWang et al.l (|2006l ), wo 
alter only the critical mass(and the k parameter simultane- 
ously) when constructing the new model, and keep the rest 
parameters describing the power law slopes and the relation 
scatter the same as those for the local Universe of SDSS. In 
this case, we now have four free parameters that are needed 
to be constrained. 



By fitting both the stellar mass function and the cor- 
relation functions of VVDS observation, we get our best-fit 



model . Stellar mass function is provided bv IPozzetti et aD 
|2007D . in the redshift bin of [0.7,0.9], with a mean redshift 
of 0.8 1. The correlation functions are from iMeneux et al.l 
(2008), based on the galaxy sample in the redshift range 
2 = [0.5, 1.2], with mean redshift of 0.85. The best fit is 
defined as the one that makes the S minimum. 

= _ x 2 W 



Xcorr 
■*V corr 



with 



x 2 m 



and 



/(corr 



N C orr 



Jog<t> — log^VVDS -, 

<r(log §vvds) 



dogw(r p ) - log w(r p ) yyps ,2 
cr{\ogw(r p ) VV Ds) 



N$ = 7, is the number of points over which the stellar mass 
function is measured, ranging from 10 9 ' 5 M© to W 116 M Q . 
N CO rr = 18, is the number of points over which the correla- 
tion function is measured, ranging from 0.2 to 10.0/i _1 Mpc, 
in three different stellar mass bins. 

Our best-fit model has the parameters: Mo — 6.31 x 
1O 11 /i" 1 M , logfc = 10.48 for central galaxies and M = 
5.03 x lO^h^M®, logfc = 9.99 for satellite galaxies. The re- 
sulting 3 = 3.07591, with x 2 ($)/A* = 1.944. These param- 
eter values are listed in Tab. 1 to be compared w ith the best- 
fit pa rameters of modelling SDSS observation (I Wang et al.l 
120061 ). Fig. [3] shows the best-fit model results. Symbols with 
error bars are the VVDS observation, and dashed lines are 
the model results. The stellar mass function is well fitted. 
The clustering is also reasonably reproduced. In Fig. [4] 
we plot the derived best-fit M a tars-M in f a u relation at red- 
shift of 0.83, for central(black solid line) and satellite(black 
dashed line) galaxies respectively. In comparison, we over- 
plotted the relations at redshift in red lines(solid line for 
central galaxies and dashed line for satellite galaxies). The 
result shows that for a given infall mass of hosting halo, the 
galaxy mass changes with redshift, and the dependence on 
redshift depends on mass. For massive haloes, the central 
galaxies inside these haloes are a bit more massive than the 
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Figure 4. The best-fit M s tars-M in f a n relation at z = 0.83 
of central galaxies(black solid line) and satellite galaxies (black 
dashed line), by fitting both stellar mass function and correla- 
tion functions of VVDS observation. For comparison, relations of 
central(red solid li ne) and satellite galaxies(red dashed line) from 
IWang et all l|2006l l at z = are also plotted. 



galaxies within the same mass of haloes at z = 0. For less 
massive haloes, the mass of central galaxies is smaller at 
higher redshift. For satellite galaxies, however, the mass of 
galaxies is much smaller toward higher redshift at all mass 

scales. 

In a recent paper, iMoster et al.l (|2009l ) claim that stel- 
lar mass function alone is enough to constrain the relation 
between galaxy stellar mass and its hosting halo mass, since 
the fit to correlation functions has a much wider range at 
its x 2 minimum than the fit to stellar mass functions. How- 
ever, this may not be true for higher redshift case. In Fig. [S] 
we show the derived stellar mass function and correlation 
functions for the best-fit model when fitting only the stel- 
lar mass function of VVDS observation. The results show 
that the clustering of galaxies is generally over-predicted 
in this case. Therefore, we believe that fitting both stellar 
mass function and the correlation functions simultaneously 
is required to get reasonable fit, at least for the current ob- 
servational data we can get. 



3 MASS INCREASE FROM Z~ 0.8 TO Z= 

From high redshift to the present day, dark matter haloes 
get larger through mergers, and galaxies inside them also be- 
come bigger in size. The galaxies gain their masses through 
either mergers with other galaxies, or by forming new stars. 
Using the model we build in Sec. [5] we already know the 
galaxy masses at redshift of 0.83. We also have galaxy stel- 
lar masses of today acc ording to the model built at z = 
from I Wang et all ((20061 ) . The galaxy mass of today is a to- 
tal amount of stellar component of galaxy stellar mass that 
already exists at redshift of 0.83, the mass increase result- 
ing from mergers with other galaxies, and in addition the 
newly formed stars during the time interval. By tracing the 
merger histories of haloes/subhaloes and hence the galaxies 



that reside in these haloes/subhaloes, the amount of stellar 
mass that was added through mergers can be calculated. 
Combined with the stellar mass of galaxies at both z = 0.83 
and 2 = 0, the stars that were newly formed during the time 
interval between these two redshift epochs can be predicted. 

For a galaxy that resides in a halo of given mass at 
2 = 0, we trace back through merger trees to its most mas- 
sive progenitor at 2 — 0.83. We plot in Fig. |S]the median 
relation between the stellar mass of the most massive pro- 
genitor at 2 = 0.83 of a galaxy and the mass of its host 
halo at present day in black solid line. Among the galaxies 
that merge into this most massive progenitor, some of them 
are galaxies that already exist at 2 = 0.83, including both 
central and satellite galaxies at that time. The other galax- 
ies are newly formed galaxies after z = 0.83, and merge 
into the main group before the present day. From our fit- 
ted model results we know that the Mstara-Minfaii relation 
evolves with time. Therefore, at the time of each merger, 
the mass of the merged galaxy should not be the same as its 
mass at the time of redshift 0.83. We get the galaxy mass 
at the time of each merger by interpolating M s t ar s-Mmfall 
relation between 2 = 0.83 and z = 0, assuming that the 
model parameters evolve linearly with redshift. 

In Fig. [6] the dotted black line is the sum of the stel- 
lar masses at 2 = 0.83 of the most massive progenitor and 
of its satellites merged in. The Red solid line is the result 
when the merged mass from central galaxies is added, in- 
cluding both the stellar mass existing at z = 0.83 and those 
newly formed since then. The contribution from the merged 
central galaxies is much smaller compared with the mass 
from merged satellite galaxies. Blue lines are the mass of 
galaxies of present day, according to our model fit result for 
SDSS observation l|Wang et al.ll2006l ). In the bottom panel 
of Fig. [5] the corresponding mass ratio of each component 
to the galaxy of the present day is shown. 

From Fig. [S] we can tell that for galaxies that reside 
in haloes of mass less than 10 12 /i -1 Mq, the mergers since 
z — 0.83 contributes to the stellar mass growth by a very 
small fraction, which can even be ignored. Compared with 
the galaxy mass at present day, their most massive progen- 
itors contribute from about 20 percent to around 60 per- 
cent of the present day mass, while the rest of the mass 
of z = galaxies should come from star formation of the 
central galaxy itself. However, for high mass galaxies whose 
hosting halo masses are more than 1O 12 ' 5 /i _1 M0, the story 
is totally different. The mass of the most massive progen- 
itor galaxy at z — 0.83 is comparable to its present day 
mass. When taking into account the stellar component of 
other merged galaxies, the total mass is significantly larger 
than the galaxy mass of the present day. This paradoxical 
result is the consequence of hierarchical merging in the cur- 
rent cosmological m odel. Here we have adopte d the merger 
trees constructed bv lde Lucia fc BlaizotJ (120071 ) for galaxies 
by taking into account the dynamical time scales. We found 
that the merged fraction of stell ar mass does not change 
when the merger time scale of ljiang et al.l (|2008l ) is adopted. 
There are several possibilities to reconcile the observations 
at low and high redshifts in the hierarchical model. One is 
that satellite galaxies are tidally disrupted in a significant 
amo unt of stellar mas s before merged into the central galax- 
ies (|Yang et alj|2009l ; IWetzel fc White! 120091 ). However, as 
we see in §4, the significant tidal disruption is not strongly 
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Table 1. Best-fit parameter values for the relations between Mi n j a u and M s t ar s in different models to fit SDSS and VVDS observations. 
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Figure 5. The best-fit model result when fitting only stellar mass function from VVDS observationl lPozzetti et al.ll2007Tl . Symbols are 
observation results, and dashed lines arc the best-fit model results. 



required by current observational data, because the unified 
model (details in §4) in which we assume no tidal disruption 
matches the observational data equally well. Another possi- 
bility is that central galaxy mass observationally determined 
may differ from the galaxy mass defined in simulation, due 
to the limit size of the observed central region of a galaxy. 
Besides, when counting for the merged mass in simulation, 
all stellar masses of the merged galaxies is added, while ob- 
servationally, the more extended stellar halo around the cen- 
tral galaxy may not be fully counted in. Therefore observa- 
tionally determined galaxy mass can be smaller than the 
simulated total mass. With the model advocated here, we 
expect to study and discriminate these possibilities when 
the model parameters can be determined better with fu- 
ture high redshift samples of galaxies. Future observations 
of intra-cluster (group) stars can also help testing these pos- 
sible mechanisms. 

In any case, qualitatively we should be able to conclude 
that for low mass galaxies, they gain quite a fraction of their 
present day mass through star formation from z = 0.83 to 
today, while for high mass galaxies, most of their present 
day mass already exist at redshift of 0.83, and star forma- 
tion is not ac tive for thes e galax ies. This is consistent with 
the result of IWang et al.l (|2007l ) where they constrain the 
star formation histories of galaxies by fitting the spectral 
distribution properties of galaxies. This is also consistent 
with the well-known fact that low mass galaxies are in gen- 
eral more active and blue in colour, while high mass galaxies 
are red and passively evolved. 

In I Zheng etafl <|2007l ). they use the HOD method to 



model the luminosity-dependent projected correlation func- 
tion of DEEP2 and SDSS surveys, and study the stellar 
mass growth of galaxies, by estimating the stellar mass from 
galaxy luminosity and colour based on the mean relation be- 
tween these properties. Compared with their result shown 
in their Fig. 9, the contribution of high redshift galaxies to 
galaxy mass of the present day has a similar trend with halo 
mass, but the absolute values of our model are in general 
higher than their result. Besides, the increase at halo mass 
of around 10 12 /i _1 Af Q is more steeper in our model. No- 
tice that they are modelling observation of DEEP2, whose 
redshift range is around 1, while our model is to fit the 
VVDS, with redshift of around 0.8, this discrepancy is in 
the right direction to be explained by t he increase of galax y 
mass through time. As pointed out bv lZheng et al. 
around 25 percent more of the stellar mass could have been 
in place in the z ~ 1 progenitors due to the fact that the 
DEEP2 sample they used are not entirely volume limited for 
red galaxies. This would decrease the discrepancy between 
these two model results. However, for massive haloes, galax- 
ies inside them are still more massive at higher redshift in 
our model than in their model result. This may be caused 
by the fact that the DEEP2 sample could have missed a 
quite fraction of red massive galaxies because of their colour 
selection of target galaxies. 



4 UNIFIED MODEL 

When a central galaxy falls into a larger group and becomes 
a satellite, the gas around the galaxy is shock-heated. The 
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Figure 7. The best-fit unified model result s of stellar mass functions and cor relation functions at z = 0(upper panel s) and z = 0,83 (lower 
panels), compared w ith results from SDSS llLi fc White! [20091 ; iLi et alj|2006l) (upper panels) and VVDS observation llPozzetti et al.ll2007l ; 
M eneux et al.ll2008h flower panels). Symbols with error bars are from observation, and dashed lines are from our best-fit unified model. 
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Figure 8. The M a tars-Mi n f a u relation for the best-fit unified model at z = 0(left panel, black lines) and at z = 0.83(middle panel, 
blue lines). Solid lines are the relation for centr al galaxies and das hed lines are results for satellite galaxies. In the left panel, the Msturs- 
Mi n fa,u relations from the model at z = of IWang et al.l l|2006l ) are over-plotted in solid green lines for central galaxies and dashed 
green lines for satellite galaxies. In the middle panel, solid and dashed green lines are relations for central and satellite galaxies from 
our best-fit model in Sec. 2 when fitting VVDS observations only. The right panel compares the median relations for central galaxies at 
z = 0.83(blue line) and at z = 0(black line). 
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Figure 6. Galaxy mass increase from z = 0.83 to z 0. for 
different hosting halo masses. Upper panel: for a given mass of 
haloes at present day, black solid line shows the median value of 
the most massive progenitors of the central galaxies that reside 
in them. Dotted black line is the total amount of mass of both 
the most massive progenitor and the merged galaxies which were 
satellites at z = 0.83. Red solid line represents the total mass of 
the most massive progenitor at z ~ 0.8 and all merged galaxies 
till z = 0. Blue line is the central galaxy mass of present day, 
which represents M s t ars -M in ^ aii relation at z = 0. Bottom panel: 
corresponding ratio of each mass to the galaxy mass of present 
day. 



gas stops cooling and the star formation ceases in a short 
time scale. The stellar mass of the galaxy can increase by 
only a small amount through star formation. On the other 
hand, the galaxy stellar mass may decrease a bit due to the 
tidal stripping effect. We now assume that in total, the stel- 
lar components of satellite galaxies do not change compared 
with the mass at the time when they are central objects. In 
this case, M s tars-Mi n f a ii relation for a satellite galaxy is the 
same as the M 3 tars-Mi n f a ii relation for centrals at an earlier 
epoch, when the satellite galaxy falls into a larger group. 

We now try to build a unified model to describe the evo- 
lution of M 3 tars-Mi n faii relation with redshift, for both cen- 
tral and satellite galaxies. We assume that at all redshifts, 
the median of the M s t a rs-Minf a ii relation can be described 
by a double power law form. We assume that the power law 
indexes and the scatter around the median relation are fixed 
with time, while the critical mass and the corresponding nor- 
malization parameter evolve with time linearly. At any given 



redshift, the stellar mass of central galaxy is connected to its 
host halo mass according to the relation. For satellite galaxy, 
its stellar mass is connected to its halo mass at infall time, 
according to the Mstars-M in faii relation at the time when 
the galaxy falls into a larger group and becomes a satellite. 

Based on this picture, we can get the stellar masses 
for all galaxies at any given redshift, and calculate their 
stellar mass functions and correlation functions. To get the 
best fit parameters that describe the evolution relation, we 
fit at the sam e time both the observation from SDSS at 
local Universe dLi fc White! |2009| ; iLi et al.l l200fj|). and the 
VVDS results at red shift of around 0.8 l|Pozzetti et"ai1l2007l ; 
iMeneux et al. 2008). For simplicity, we fix the slopes de- 
scribing the M s tars-Mi n f a u relation to be the same as the 
slopes of t he relation of cen tral galaxies when fitting SDSS 
data onlv l|Wang et al.ll2006l ), i.e., a = 0.29 and (3 = 2.42. 
In total, we now have 5 parameters to fully describe the 
M s tars-M in f a u relation: critical mass Mo and ko for red- 
shift 0, Mo. 8, ko.s to describe the relation at z = 0.83, and 
the scatter a around the median relation. 

The best-fit model is determined when the resulting 3 
gets its minimum. 
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N<s>,vvds = 7, N<s> : sdss = 29 are the numbers of points 
over which the stellar mass functions are measured in two 
observations. N corr ,vvDS = 6x3, N CO rr,SDSs = 30 x 5, are 
the numbers of points over which the correlation functions 
are measured in different stellar mass bins. 

Our best-fit model has the parameters: Mo — 3.21 x 
1O 11 /i _1 M , logfco = 10.17, Mo.s = 4.34 x lO^ft^M®, 
logfco.s = 10.15, and a = 0.240. The resulting 3 = 14.455, 
with X 2 {®o)/N$, S DSS = 3.133, xlorr,o/N C orr,SDSS = 3.218, 

X 2 {®0.8) /N$,VVDS = 3.778, and xlorrfi.&l Ncorr,VV DS = 

4.325. The parameters are also listed in Table. 1. Fig.[7]shows 
the best-fit model results, compared with the observation 
from SDSS and VVDS surveys. All the observations are rea- 
sonably reproduced, although the resulting \ 2 is larger, and 
the fit is a bit poorer than the previous model in Sec. 2. The 
constrained evolution of critical mass M and normalization 
parameter k can be described as: 



M(t) 



*(*) = 



t -t 

to — to. 8 



x(M .8 - M ) + M 



t -t 



x(k .8 - ko) + ko 



to — to. 8 

t is the age of the universe at different redshifts. to = 
13.6098G?/r, and to. 8 = 6.7531Gyr, are the age of the uni- 
verse at redshift and 0.83 respectively. 

Fig- El shows the best fit Mstars-M in f a u relation at two 
different redshifts. For comparison, we also over-plot the 
model result when fitting SDSS and VVDS data separately 
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in green lines. We find that at both high and low redshifts, 
central galaxies in massive haloes are less massive in the uni- 
fied model than the separate model. The difference of satel- 
lite galaxies between the two models is smaller. Although 
the best-fit M s tars-Mi n f a ii relations are different in our two 
models, they can both give reasonable fits within the obser- 
vational error range. The main difference of the two mod- 
els exists for massive galaxies. Considering that the number 
of high mass galaxies is small in observation, which causes 
large error bars on the data points, especially for VVDS 
survey, the statistics of high mass galaxies are not tightly 
constrained. Therefore different sets of parameters can both 
give reasonable fits to the observational data. The models 
can be better constrained only with improved observational 
data. 

In this best-fit unified model, from Fig. |H]we know that 
both at z = and z = 0.83, satellite and central galax- 
ies have similar relation for massive haloes, while for low 
mass haloes, central galaxies are more massive than satel- 
lites. In the most right panel of Fig. [8j we compare the 
M s tar s-M i n f a u relations for central galaxies at redshift 
and 0.83, for the best-fit unified model results. The evolu- 
tion of this relation is larger for galaxies in low mass haloes 
than in more massive haloes. For all masses of haloes, cen- 
tral galaxies that reside in them are more massive at lower 
redshift than at higher redshift. 

Based on the method studying the mass increase of 
galaxies from redshift of 0.83 to today described in Sec. [31 
we can also detect the mass increase situation according to 
the result from the unified model stated above. We perform 
the same analysis as previously stated, and show the re- 
sults in Fig. [9] We find that the most massive progenitors at 
z = 0.83 contribute from around 25 percent at halo mass of 
~ 10 11 hT 1 Mq to around 70 percent for haloes massive than 
~ 10 12 ' 5 /i _1 Mq to the galaxy mass of present day. This frac- 
tion is much smaller for massive galaxies than the separate 
model. The merged fraction of galaxies, however, is higher 
in the unified model. Nevertheless, qualitatively the same 
as in the separate model, the unified model also shows that 
for galaxies in low mass haloes, they gain their mass from 
z = 0.83 to today mainly through star formation, and extra 
mechanisms are needed to explain the excess of the total pro- 
genitor masses to the present-day mass for massive galaxies. 
Since in the unified model the tidal tripping is assumed to 
be unimportant, the excess of the stellar mass is presum- 
ably located in the stellar haloes of central galaxies. With 
upcoming larger deep samples of galaxies, we will consider 
to include the tidal stripping effect in the HOD modeling in 
a future work. 



5 STELLAR MASS FUNCTIONS AT HIGHER 
REDSHIFTS 

In Sec. [2] and Sec. [4] we build models based on the observa- 
tions of SDSS and VVDS data of galaxy samples of the local 
Universe and of redshift of ~ 0.8. Assuming that our models 
can be applied to even earlier epochs, with the parameters 
evolving linearly with redshift, we can predict stellar mass 
functions of galaxies at higher redshifts. 

Fig. [10] presents the results of stellar mass functions 
at different redshift bins. Black solid lines are results from 
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Figure 9. The same as Fig.[6]showing galaxy mass increase from 
z = 0.83 to 2 = for different host halo masses, except that for 
the unified model results. 



recent observation of iKaiisawa et al.l (|2009h , for four red- 
shift bins. The stella r masses are normalized to the Chabrier 
IMF (|Chabrierll2003l '). The median redshift of 0.5 < z < 1.0, 
1.0 < z < 1.5, 1.5 < z < 2.5 and 2.5 < z < 3.5 
galaxy samples are 0.802, 1.172, 2.080, and 2.952. In the 
upper left panel, also plotted is the result of lPozzetti et al.l 
l|2007t ) of VVDS observation in redshift range of (0.7,0.9), 
which is used in the previous sections to constrain our 
best-fit models. I n the l ower two panels, observations from 
iMarchesini et all (|2008l ) are over-plotted, for 1.3 < z < 
2.0(green symbols) and 2.0 < z < 3.0(red symbols) galaxy 
samples. Blue dotted and dashed lines are predictions from 
best-fit models presented in Sec.[2]and Sec.[4j at redshifts of 
0.83, 1.17, 2.07, and 3.06. The results from different obser- 
vations are consistent with each other, and the predictions 
from both of our models are also consistent with these ob- 
servations. 

Notice that at all redshifts, the prediction from the uni- 
fied model gives a higher amplitude of stellar mass function 
at low mass end, and a lower amplitude at high mass end, 
compared with the prediction from the model where SDSS 
and VVDS data are fitted separately. In the unified model, 
galaxies have smaller masses at high mass end of haloes, and 
the satellite galaxy mass is higher at high redshift than the 
separate model. These together cause the difference of am- 
plitudes of stellar mass functions at low and high mass end. 
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As discussed in Sec. 4, although the predictions from our two 
models differ at all redshift bins due to the obvious differ- 
ence of the M star s-M 'i n fail relations obtained in the models, 
it is hard to tell which model is a better description of the 
real universe observed, due to the poor statistics of current 
data. Observations with smaller errors will be able to help 
us build a more precise model. 



6 CONCLUSIONS AND DISCUSSIONS 



We apply the empirical model of IWang et al.l (|2006T ) to 
higher redshift, to link galaxy stellar mass with its host- 
ing dark matter halo mass at redshift of around 0.8. The 
M s tars-M in f a u relation is constrained by fitting both the 
stellar mass function and the correlation functions at differ- 
ent stellar mass intervals from VVDS observation. We find 
that the difference of M s tars-M in f a u relation between cen- 
tral and satellite galaxies at z ~ 0.8 are much larger than the 
galaxies at present day. At all mass scales of haloes, satel- 
lite galaxies at z = 0.8 have much lower stellar masses than 
satellites within the same mass haloes at z = 0, while the 
mass of central galaxies changes less from high redshift to 
today. Central galaxies are less massive in low mass haloes, 
and more massive in high mass haloes at z = 0.8 than the 
galaxies at z = 0. 

Under the assumption that the mass of satellite galaxy 
remains about the same as it is a central before it infalls to a 
larger group, we build a unified model for the M s tars-Mi n f a ii 
relation, which describes the evolution of galaxy mass as a 
function of halo mass at any given redshift. Satellite stellar 
mass is determined by the M s tars-Mi n f a u relation of cen- 
tral galaxy at the time of its infall. The best-fit model of 
both SDSS and VVDS stellar mass functions and cluster- 
ing functions gives an obvious evolution of M gta rs-Mi n faii 
relation from z = 0.8 to z = 0, with the mass of galax- 
ies at higher redshift being lower than the galaxy mass at 
the present day, for all masses of hosting haloes. The M s t ars - 
Mi n f a ii relation provided by the unified model is different as 
the relations in the models when SDSS and VVDS data are 
fitted separately, especially for galaxies with massive haloes 
and at high redshift. The central galaxies are less massive 
and satellite galaxies are more massive in the unified model 
than the separate model. Different sets of parameters can 
both give reasonable fits to the observed data, because the 
statistics of high mass galaxies are not tightly constrained 
in observation at high redshift. 

We study the amount of galaxy stellar mass growth 
from 2 ~ 0.8 to today, in either way of galaxy merger or 
star formation. For the models when SDSS and VVDS data 
are fitted separately, we find that for galaxies that reside 
in haloes of mass less than 10 12 ft _1 M@, the masses of their 
most massive progenitors at z — 0.83 vary from about 20 
percent to 60 percent of the present day mass. Meanwhile, 
galaxy mergers contribute only a small fraction to the galaxy 
mass of today. This indicates that a large fraction of z = 
stellar masses comes from star formation during the period 
between z — 0.83 and z = 0. For galaxies within massive 
haloes, the total amount of stellar mass from the main pro- 
genitor at z = 0.83 and from that increased from mergers 
with other galaxies actually exceeds the present-day galaxy 
mass. Although there could be a difference of the definitions 



for stellar mass of central galaxies in observation and in the 
models based on simulation (for example, the stars in the 
envelope of the central galaxies may not be counted in obser- 
vation), this indicates that there may be little star formation 
ongoing in these massive galaxies. The unified model basi- 
cally tells the same story, except that the contribution of 
the most massive progenitors at z = 0.83 to the mass of the 
present day is no more than 75 percent even for the most 
massive galaxies. 

Based on the models we build, we give predictions of 
higher redshift stellar mass functions, with redshift up to 
z ~ 3. It is encouraging that our predictions are consistent 
wi th recent observations o f ste llar mass functions pre sented 
bv lMarchesini et all (|2008l ) and lKaiisawa et al.1 (120091 ). How- 
ever, the amplitude of the stellar mass functions predicted 
by the unified model is always higher at low mass end, and 
lower at high mass end, compared with the prediction of the 
model that fits SDSS and VVDS data separately. Besides, in 
both of our models, we assume that for the M s tars-Mi n faii 
relation, its slopes at high and low mass end and scatter 
around the median value do not change with time, which 
may not actually be true. The current models can be tested 
by future observational data, and a more accurate model 
can be constructed with better data. 



Besides the observational data of VVDS survey used 
to constrain our model in this work, galaxy stellar 
mass function and correlation functions in different stel- 
lar mass bins are al so o btained for zCOSMO S survey by 
iPozzetti etafl d2009h andlMeneux et all (|2009h . As pointed 
out in lMeneux et al.l (|2009h . the zCOSMOS field is centred 
on an extreme overdensity region. We therefore choose to 
use the VVDS results in this modelling work, although it 
is possible that the VVDS field is a bit under-dense on the 
other hand. Hopefully this cosmic variance effect can be sig- 
nificantly reduced with upcoming larger samples of galaxies, 
based on which the models can be constrained much more 
tightly. 



Recent study of iMoster et al.l l|2009h determined the 
galaxy stellar mass - halo mass connection at high redshift 
constrained with stellar mass functions only. They claim 
that stellar mass function alone is enough to constrain the 
relation between galaxy stellar mass and their hosting halo 
mass. We have shown at the end of Sec. [2] that this may 
not be true for high redshift situation, although this was 
proved viable when considering models at the local Uni- 
verse. However, their result shows the same trend as our 
M s tars-M in f a u relation in the separate model, although the 
quantitatively value of the relation is somewhat different. 
For low mass haloes, galaxies at higher redshifts have less 
stellar mass than galaxies at a lower redshift. For high mass 
haloes, galaxies at high redshift is a bit more massive than 
galaxies at low redshift. In our unified model, on the other 
hand, galaxy mass is always smaller at high redshift than 
that at low redshift, but the difference for galaxies within 
massive haloes is quite small. The precision of all these mod- 
els can be tested only by improved observational results both 
on stellar mass functions and on the clustering properties at 
high redshifts. 
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Figu re 10. The pr ediction of stellar mass functions at high redshift of our best-fit models. Black points are observations from 
iKaiisawa et al.l 1120091) , for four redshift bins. Blue dashed lines are the predictions when applying the M s tars-M in f a u relation of our uni- 
fied model to higher redshifts of 0.83, 1.17, 2.07, and 3.06. Blue dotted lines are predictions of the best-fit models when S PSS and VVDS 
data a re fitted separately, assuming parameters evolve linearly with redshift. In the upper left panel, the observation of Pozzctti ct al. 
(2007) is plotted(blue sy mbols), for 0.7 < z < 0.9 galaxy sample, which is used for building our model at z0.83 in Sec. [2] In the lower two 
panels, observations from Marchesini et al. | l|2008h are plotted, for 1.3 < z < 2.0( green symbols) and 2.0 < z < 3.0(red symbols) galaxy 
samples. Stellar masses of galaxies are normalized to the Chabrier IMFljChabricr 20031). 
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